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When the equatorial spin velocity, ti, of a charged conducting sphere approaches c, the Lorentz 
force causes a remarkable rearrangement of the total charge q. 

Charge of that sign is confined to a narrow equatorial belt at latitudes b ^ \/3 (l — v^/(?) ' while 
charge of the opposite sign occupies most of the sphere's surface. The change in field structure is 
shown to be a growing contribution of the 'magic' electromagnetic field of the charged Kerr-Newman 
black hole with Newton's G set to zero. The total charge within the narrow equatorial belt grows as 

(l — v^/(?^ ^ and tends to infinity as v approaches c. The electromagnetic field, Poynting vector, 
field angular momentum and field energy are calculated for these configurations. 

Gyromagnetic ratio, g-factor and electromagnetic mass are illustrated in terms of a 19th Century 
electron model. Classical models with no spin had the small classical electron radius /mc? ~ 
a hundredth of the Compton wavelength, but models with spin take that larger size but are so 
relativistically concentrated to the equator that most of their mass is electromagnetic. 

The method of images at inverse points of the sphere is shown to extend to charges at points with 
imaginary co-ordinates. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

When a charged conductor rotates steadily, any small 
resistivity in the conductor will cause dissipation, un- 
less the charge rotates with the conductor. Thus, in the 
steady state the surface current is due to the rotation of 
the surface charge. For a sphere rotating with equatorial 
speed V much less than c, the charge is almost uniformly 
distributed over the surface and the resulting magnetic 
field is almost uniform within the sphere and externally 
dipolar. If the charge were frozen to be uniform over 
the surface and the sphere were an insulator then those 
would be the forms of the fields, however fast the sphere 
rotated. (Appendix 1). 

By contrast, when the sphere has a conducting surface 
the Lorentz force pushes the charge along the surface to- 
wards the equator, but the effect is only pronounced at 
relativistic speeds. Beyond v = 0.93c the effect is so 
strong that the charge density at the poles is of the op- 
posite sign to the imposed total charge q. As v gets 
still nearer to c more and more of the sphere is cov- 
ered by such reversed charge, while the charge of the 
same sign as q is confined to an ever narrowing belt 
around the equator. However this belt now contains a 
charge considerable greater than q since the excess is 
needed to cancel the reversed charge elsewhere. Indeed 

we show that as v approaches c the total charge in the 
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belt becomes ^ (l — v^/c^) * which tends to infinity. 
The bound aries of th e equatorial belt are at latitudes 
|5| ~ \/3 -^/l — v'^/c^ so the width of the belt tends to 
zero. 

Defining u by (1 — u)/{l + u) = \/l ~ v^jt? we find 
that as V increases, the external field is changed by the 
addition of the electromagnetic field of a charge at the 



imaginary point (0, 0, w-^/u). This is the Magic electro- 
magnetic field discussed earlier (Newman, 1973, 2002; 
Lynden-Bell, 2003) which may also be found as the field 
left by a charged Kerr-Newman black hole (Newman et 
ai, 1965) when G is set equal to zero. In our context 
it is being generated quite naturally as w — > c. The 
corresponding internal field is that due to a charge at 
the complex image point under inversion in the sphere 
(0, 0, ~ia/y/u). Wave equations for particle motion in the 
Magic field have remarkable separability (Carter, 1968b; 
Chandrasekhar, 1976; Page, 1976; Lynden-BeU, 2000). 
All the Kerr-Newman metrics have the same gyromag- 
netic ratio as the Dirac electron (Carter, 1968a) and so do 
all the conformastationary metrics and others remarked 
upon by Pfister and King (2002). 

There are several points of principle that arise in this 
problem which we raise but do not discuss in detail. 

1. If a real physical sphere with a conducting sur- 
face is rotated relativistically, what shape would 
it become? One might suppose that it ought to 
be a sphere in the rotating axes, but the metric 
in such axes is anisotropic, so it is not clear what 
shape should be called spherical in such axes. We 
sidestep all such discussion by taking the surface 
to be spherical in fixed axes however fast the body 
is rotated. By this we elimiate all discussion of 
Fitzgerald contraction, the physical constitution of 
the sphere and its distortion under centrifugal and 
Lorentz forces. While this may reduce the physi- 
cal reality of the problem, it gives a well defined 
uncomplicated problem. 

2. For a slowly rotating conductor the charge resides 
on the surface so it does not matter whether the 
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whole body of the sphere is a conductor or only a 
surface layer. However, once the magnetic part of 
the Lorentz force is significant, the condition that 
the Lorentz force density must vanish inside the 
sphere leads to a requirement for internal charges 
and the currents due to their rotation. Thus in 
the relativistic regime the conducting solid sphere 
gives a problem that differs from the conducting 
spherical shell problem solved here. 

3. We do not assume that our conductor has zero resis- 
tivity, indeed we rely on some resistivity to ensure 
that in the steady state the surface charge rotates 
with the body. In particular, we do not assume 
that either the surface or the body of the sphere is 
a superconductor, so questions of the expulsion of 
magnetic flux by the Meissner effect do not occur 
in our problem. 

4. A number of very interesting problems arise if it 
is assumed that the charge carriers in the conduc- 
tor have mass, because at relativistic speeds the 
inertia of the charge carriers will itself cause them 
to concentrate toward the equator. However the 
Maxwell equations themseves already account for 
that part of the inertia contained in the electro- 
magnetic field. When an electron is accelerated 
the electromagnetic field is accelerated with it, so 
at least part of its mass is due to the inertia of its 
field energy. To put in the full mass of the electron 
as well as the Maxwell stresses of the field distri- 
bution would involve some double counting, not to 
mention the larger effect of inertia on the heavy 
nuclei that form the conductor's lattice. We shall 
again sidestep these problems by assuming that the 
Lorentz force component along the conducting sur- 
face is zero, so that there is no effect of inertial 
forces on the charge carriers (other than those on 
their electromagnetic fields as implied by Maxwell's 
equations). 



II. THE MAXWELL EQUATIONS FOR THE 
PROBLEM 



We call the radius of the sphere a and rotate it with 
angular velocity Vl so Via = v. We shall find it convenient 
to use ljo = fla/c = v/c. In the steady state, the Maxwell 
equations reduce to Curl E = 0, divE — Anp, divB — 
and CurlB = inj. On the sphere's surface there is a sur- 
face charge a{9) and a sheet current = anac~^ sm6 
but away from that surface there are neither charges nor 
currents so we may write E = — V<i> and B = — VX where 
V^<1> = = V^X. The z —z symmetry shows E even 
and B odd. 

The solutions to Laplace's equation obeying the 
boundary conditions at r = and oo are with /i = cos d 



/Qs 2n-|-l 



r 

J- \ 2n 

a 



and 



1 /a\2«+2 



^ 2n + 2 



y-j ^2n+l^'2n+l (P-) T > a 



a sr 



1 /r\2n+l 

1 \a) ^27i+iP2n+i (m) r < a 



^ 2r 



Here q is the total charge on the sphere and the co- 
efficients have been arranged so that $ and dxjdr are 
both continuous across r = a. This corresponds to the 
requirements that the surface component of E and the 
normal component of B must be continuous. 

From the discontinuity in Er we have the expression 
for the surface density of charge 



Ana = qa ^ ^ (4n 1) $2n ^2« (a^) 



(1) 



so the condition that q is the total charge gives $o = 1 
on integration over the sphere. The discontinuity in Bg 
gives the surface current so 



47rj0 = 

E (2»+2K2Ui) ^2»+iP:^„+i (/^) , (2) 

where Pan+i = d/d^i{P2n+i). 

The condition that the current is solely due to the ro- 
tation of the surface charge combines the above equations 
to give 



E 

n=0 



in + 3 



{2n + 2){2n+ 1 



-^2n+l P2n+1 



OJ 



^ (4n + 1) $2« P2n 



(3) 



Our final physical requirement is that at equilibrium 
there is no surface component of the Lorentz force (be- 
cause the surface is a conductor) 

Eg+u! sin 6 Br = , 
which implies in terms of our potentials 



oo 



oo 



^ ^2«-P2n M = E ^2n+lP2n+l (m) ■ (4) 



Our problem has now been reduced to the mathemati- 
cal one of solving equations ||2Jl and Q for the coefficients 
^2n and X2n+i when w ^ 1 and is given, and $o = 1- 
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III. MATHEMATICAL SOLUTION 

We need the following mathematical properties of the 
Legendre Polynomials 



might neglect the 6m for m ^ 2. This will form the basis 
for a rapidly convergent approximation scheme. Using 
to eliminate the $ from ((HJ we find for m ^ 1 



P-ln ^2m+l 



m ^ n 
m > n 



m ^ n — 1 
m > n — 1 



together with their orthogonality relation. Multiplying 
Q by ^P2m+i and integrating from —1 to +1 we find 
using the above 



$2n = -^?72m+l 



(5) 



n— m+l 



where we have defined '72m+i as X2m+i/(4TO + 3) because 
it proves to be a more convenient variable to use. Simi- 
larly multiplying © by ■^P2m and integrating we find on 
writing the result in terms of 77 



(1 + Sm) ^2m+l 



(Jl2m+'i " 2'?2m+l + ^2m-l) 



hence we obtain a recurrence relation with almost con- 
stant coefficients 



'?2m+3 + 2 [2U-^ (1 + 5m) - 1] '?2m+l+'?2™-l = . (10) 

Thus ^h+2 [2^-2 (1 + ^) _ 1] ri^j^-q^ ^ Q. Now, since 
^ 1 , 2ti)^^ — 1 ^ 1, so with a less than 2% error in the 
middle coefficient we can set (5m = for m > 2 . H1U|) is 
then reduced to a recurrence relationship with constant 
coefficients. 
The solution is 



'?2rn+l = C(-U)" + D{~U)-"' m^2, (11) 

where u and are the roots for t of 



^ (An + 3)^ 
^ (2n + 2)(2n-|-l) 



I2n+1 



(6) 



f - 2(2w"2 - l)i + 1 = , 



notice that 



(12) 



We shall give two different methods of solving © and 
© which are complementary in that one gives enlight- 
enment on what is happening in the other. 



A. First Method 

Subtracting from 10) the same equation with m + 1 
written for m we obtain 



1 - u 



= 4u/(l + uf and ^JX-uS^ = \ 1 - ^ 

' ^ ' V 1 + u 

(13) 

Evidently the product of the roots is one and the sum 
of the roots is positive, so both roots are positive with 
u < 1, > 1, for uj < 1. However, the expression 
((nTl will diverge oscillatorily due to the term unless 
D ^ 0. Hence, 



^2m+2 = (^2m+3 " '?2m+l) , (7) 

performing the same operation on @ and dividing by 4 



{m + l){m+ ^) 



Y7''^2m+l 



-($2m-H2 - $2m) ■ (8) 



Now the coefficient of '72m+i is remarkably close to 1: 
indeed it is 



_9_ 

16 



= 1- 



16(to+ 1) (m+ i) 



= l + Sn 



(9) 



where 5q = ^, Si 



— L A^J_ A=J_ A^J_ 
g J ^ 48 ' "2 120 ' 3 224 ' "4 350 

etc, thus the Sm quite rapidly become small compared to 
one which suggests that if only 1% accuracy is needed we 



^2m+l = Ci-UY 



where, 



u = 2uj 



Notice that 



but as 



In particular 



1 - - 



m^2 



1, l-2v/l-cj2 



^5 = Cu^ , 



(14) 



(15) 



4 



and 



2 



V2n 



c- 



1 



(16) 



We still need to determine the constant C as well as 
^3, Vi and the $2ri- Our strategy is to return to the exact 
equations © and (0) for to = 0, 1, 2 and to combine them 
with the results l|15|l and l|16|l to determine the unknowns 
$2, $4, ^3, ^1 and C. 

We rewrite © 



^ (1 + S^) V2 



n+l 



(17) 



finally putting expressions (|23|l and H24() for and Vi 
into H21|) we have our equation for the constant C again 
using 1(1^ 



(25) 



where 



A = 1 + (1 + m) 6q + 5i{1 + u + u^) +5q 5i{l + uf , 

(26) 

with (^0 = I and Si = and u given by itT^ . 

With C so determined is given by H24|) and by 
()23|l may be incorporated into the more general formula 



TO = gives C- + (1 + Si) V3 + (1 + So) = uj/i , 

I + u 

(18) 



TO=1 gives C-^ + {l + Si)V3 = ^<^2 , (19) 
1 + u 4 



u UJ 

m ~ 2 gives C- = —$4 , 

^ 1 + u 4 



with TO = gives $2=^(^3-^1) 



(20) 



2m+l 



for m ^ 1 . 



For u < 0.9 a somewhat more accurate solution may 
be shown to be 



'?2m+l = C{-uy 



jl + u) 1 
16(1 -u) (to+I) 



Its accuracy is 1.5% at u = .8 but diminishes to 7% 
at M = .9. Beyond that we give a separate treatment in 
section 3.2. 

With all the known the <i>2m are given by (0 with oj 
given by (jTSJl. Thus 



with TO = 1 (71 gives ^4 ~ lo {Cv^ ~ 7/3) 
So subtracting ((T^ from ((TH|l 



$0 = 1 



$2 = -m[1 + (5i(1 + w)]/A 



(l + ^o)ryi = _[l-c^(r/3-??i)] , 



while fH^ and become 



C- 



1 + M 



and c(^^^uA^-^V, 
using ifT^ to re-express ^ as (^i^^yi we find 

putting 1(221) into we find 



= [1 + Siil + u)^] C , 



(21) 



(22) 



(23) 



(24) 



2m 



-u)™/A 



TO ^ 2 



Hence 



00 



A2 



2 r 



(1 



+ 2Si{l + u) + Sl{l + 



.(27) 



Those who wish for a much more accurate calculation 
including the first order perturbation theory for the terms 
neglected above, will find it in the appendix. However, 
there is always some tension between simple explanation 
and a detailed calculation in which the reader may get 
lost in the mass of formulae, so I have chosen the path 
of simplicity for the main paper which helps to keep the 
physics to the fore and to suppress the minutiae of math- 
ematical manipulation. The above approximate solutions 
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for the potentials are 



Now 



1 (l - ±) + ± 

-(5l?/(l+u)2^P2(/i) 



n=p ^ ' 

r> a . 



1 

aA 



n=0 ^ ' 



r < a 



(28) 



within the square bracket the first term is clearly the 
magnetic field due to the charge at the imaginary point 
i^/u az with the convention that E + iB = —V [$ + iX] 
this fits nicely with our result for the potential. The 
second term in the square bracket may be interpreted at 
the magnetic potential of the charge distribution ^ ^ 
at i^/u a'z, between a' — and a' — a. However, there is 
no corresponding part in the electrical potential. 



B. The Limit v = c 



(r2-2^6r + 62) ^ = ^ J] (^^) , r ^ \b\ . 

(29) 

So if we put b = ia^/u, the real part is 

which should be compared to the first of the sums in $. 

We see that our potential is that due to charges ^q/A 
at the imaginary points z = ± ia^/u plus a quadrupolar 
correction term in &1P2 and a charge q{l — at the ori- 
gin. Whereas that construction gives us the field outside 
the sphere, the field inside the sphere can be thought of 
as arising from charges at the imaginary points conju- 
gate to them under inversion in the sphere r = a plus a 
quadrupolar correction term in 61 P2. 

The magnetostatic potential X takes the form for r ^ a 

2n-F I /aN2"+2 



A 2n + 2 



We notice that putting b ~ i^fua and taking the imag- 
inary part of H29|l yields for r ^ a 



also 



Xm / (r^ — 2iii\/u ar 
Jo 



M2 ^ (-u)" /a\2n 

a ^ 2n + l\r) 
Hence at least for r ^ a, X is of the form 



+2 



X = 



A 



1 



— 2ifj,^/u < 
~ 2ifi\/u w 



i da 



2 Jo ' 

3(5i(1 + m)2 aui 



4A 



-q- 



(30) 



When v/c is within half of a percent of 1 (ie 7 > 20) 
even the amazingly good perturbation theory of the Ap- 
pendix which allows linearly for all the terms we have 
neglected, breaks down; it generates terms in | ln(l — u) 
which can no longer be treated as perturbations; we 
therefore need a new approach for that very small region. 
Our exact difference equation for V 2,1+1 is 



{u + lf 



l&u (?i+i)(n+l) 



, 



at u = f the solution when the term is neglected is 
'^2n+i = C(— f )". To solve the equation without any such 
neglect, we write '?2n+i = (— l)"fn and obtain 



{Vn+l - Vn) - {Vn - Un-l) - 
(u-f)2 (W+1)2 f 



16w {n+^){n+l) 



Vn^O. (31) 



Let us look for solutions in which w„ varies sufficiently 
slowly with n that we may approximate Vn as a continu- 
ous function i;(n). Then {Vn+l~Vn) — {Vn—Vn-l) — 



Thus to sufficient accuracy we have for large n 
'{u-lf {u+lf 1 



16m 



v{n) . (32) 



The solution that converges at large n is a Bessel func- 
tion of imaginary argument. Setting N = n + j 



l-u 



-N 



U2 



where 



1 (u + lf 



v{n) ^ C7r~5 r(j^) ( ^Y^^ 
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for small and to Cexp (^—^-^N^ for large. 

In the relevant regime with u ~ 1 then i/ = . 

It is more helpful to derive the correct form for n large, 
as we shall do presently. However, before that we should 
assure ourselves that our technique works by applying it 
to the equation in which the (n + |) terms are absent. 
Then the exact convergent solution is w„ = Cw" and near 
u = I this may be re-expressed as 



We have checked by computer that the asymptotic for- 
mula with these values of A and B works down to n = 
with an error with less than one part in 1000, so 

V2n+i = Ci(-l)"Af-" [1 + AN-^ + BN-^] (33) 

where = n + | is the accurate convergent solution 
with the above values of a, A, B. The constant Ci is 
determined via equation © with m = 0,<I>o = l,w = l 



Vn = Cexp{nln[l — (1 — u)]} = 
Cexp [-n{l - u)] + 0{1 - uf 

Equation (|32|l gives a solution 



w(n)=Cexp — n(l — u)u 2 



which is the same to order (1 — m)^. Thus our ap- 
proximation only works in the neighbourhood of u = 1 
but that is where we need it. When u = 1 equation 
has the exact solution ^(n) = Ci (n -I- |) " where 
a = 1/ (2\/2 -I- 2) ~ 0.207107 so this must be the asymp- 
totic form of the solution of the difference equation for 
Vn when n is large. 

Writing N = n + j the exact equation at u = 1 is 



Vn+l 



16 ) 



to be exact away from it = 1 an extra term 



1 + Ti(]^ 



Vn is needed on the right 



We look for solutions of the form 



Vn^CiN-"{l + AN-^ + BN-'^) . 

Terms of order A^^^^" in the recurrence relation give 
us the value of a found above and the terms of order 
j^-i-a ^^-^ ^gg(j gjyg equation for A. This is 




(-1)"7V-" 



[\ + AN-^ + BN-'^] 



0.699 



.358 



(34) 



The convergence of the sum is slow 2x10^ terms gave 
.6926 but comparison with 2 X (10\ 10^, 10^ 10^ 
10^, 10^, 10*) terms show the increments decrease in ge- 
ometrical progression by a factor 1.608, so extrapolation 
is accurate. 

In this regime we can now calculate Vi = .375 = |. 



IV. EXPLICIT EXPRESSIONS FOR 
POTENTIALS, CHARGES AND CURRENTS 

Using (|29|l we write expression (I28|) for $ in the form 



r ^ a 



l--(l + i^T^2(/i) 



^i?e (a^ + 2iu + 2 arjj, — ur"^ 



, r ^ a . 

(35) 

Notice that the points ^0, 0, iu'^aj and 

^0, 0, —iu^'^aj are inverse points for the sphere. 

The surface density of charge is given in terms of the 
discontinuity in d^/dr at r = a . 



^-"96 



2(2a-|-3) 



'0.0193078 . 



We chose the form with only even powers, as then the 
odd powers cancel. The next terms are of order A''"^ 
times the leading terms. These give 



47rcr = qa " <j 1 - ^ ( 1 -I- ^m(1 -|- ujP^ ] + 



— (1 + u)Re (1 + 2m3^ - u 



Putting /i = we obtain the values on the equator 



(4a-flO)B = --L(l30a + 134^) + 



1 

2048 



A 
120 



B 



0.0068443 . 



47ra2 



1 / 5 , 

1 - — 1 u 1 

A V 96 ^ 



1 + u 



A{l-u)i 



(36) 
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this clearly diverges as u approaches 1 because A re- 
mains finite. Similarly putting /i = 1 and evaluating 
. 1 

1 + iU2 



By then most of the sphere's surface is covered with 
charge of the opposite sign leaving only the equatorial 
belt with a much enhanced charge of the same sign as q. 
The charge in that belt is approximately 



'^polc 



1 - 



l + —uil + u) ] + 



1 - 3u 



A(l +m)2 



(37) 



This density changes sign as u increases. Taking ac- 
count of the variation of A given in , we find (Xpoic — 
at M = 0.454, 7 = 2.66, uj — 0.93 . The surface current 
is solely due to the advection of the surface charge so 
J0 = <7uj sm9. 

We may get a good idea of the value of 9 or (/z) at 
which the charge density changes sign as follows; is 

not far from 1. In fact 1 — ^ = where by (|26|l 



so, varies between i and ^ . The major variation 
a with ^ is caused by 



l + u 
A 



i?e I 1 + 2m 2 /i — u 



which becomes > 6 for (1 — u) < i. We write 



Me 



i0 



where 



13 = i&ir' 2u2 



and 



(1 



l-u 



2 IT 

3 2 



Evidently M 2 cos ( will become zero when /3 = 
^ i.e. where 2u^^/{l — u) ~ tan ^ — \/3, so this cor- 
responds to /i = 5-^/^(1 ~ which becomes small when 
u approaches one. Thus the charge of the same sign as q 
becomes confined to a belt between latitudes 



q{l + u) 



2y2A(l -u): 



which tends to 00 like 2A ~ 02) these only hold 

when V is close to c. A rough approximation to the charge 
in the belt more generally is 



1 / i 



7o )U{l- 70) 



where U is one for positive values and zero otherwise and 
70 is the value of 7 at which the charge changes sign at 
the pole ie 70 = ^ = 2.66 . 

Although all our calculations thus far have been in 
terms of the magnetostatic potential X we shall find that 
several global properties of the electromagnetic field are 
more nicely expressed in terms of the vector potential A. 
Since B lies in meridional planes and is axially symmet- 
rical, we may take A in the form —Q V cf) where (j) is the 
azimuthal angle of spherical polar co-ordinates. Then 



B = Curl A = -VC X \/<f) 

Br = -^X/^r^-^^^^Q/^e^+r-^^C|^^l . 
Inserting our expansion for X in terms of the P2n+i 



q 



2a? ^ 



— ^2n+l -'2n+l — , 

r / a/i 



but the Pn obey Legendre's equation 



1 



(2n+l)(2?i + 2)d^ 



2n+l 



P 



2n+l 



SO inserting this into the expression above and using the 
boundary condition that C = when /z = -1-1 we find 



2"+l X2n+1 

rJ (2n-|-l)(2n-f 2) 



for r ^ a and a very similar condition for r ^ a. 
Now 



± sin ^ 



and the total width of the belt around the equator 
\/3a(l — u). 



and 



(1 - fi'^)dPn/dfi = n(P„_i - fiPn) , 
^iPn = [(n + l)P„+i + nP„_i]/(2n + 1) 
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V. GLOBAL PROPERTIES OF THE 
ELECTROMAGNETIC FIELD 

The Poynting vector is x B and as this flows 

around the axis it gives an angular momentum density 



— r X (E X B) . 

Aire ^ ' 

Thus the total angular momentum in the electromagnetic 
field is 



L= — /rx(Ex B)rf^a; = +— / r 

x[V$ X (VC X V(f))]d^x , 

where we have used axial symmetry and the fact that B 

lies in meridional planes to write the vector potential in 
the form — CV^ where (j) is the spherical polar coordinate. 
Since = we deduce that 



FIG. 1: The electric lines of force for a relativistically ro- 
tating conducting spherical shell with u=0.9 i.e. gamma=19, 
v/c=0.9986. Notice the charge concentration at the equator 
and the diminution of the field strength as the lines cross the 
sphere again, due to the opposing charge there. In contrast 
the static shell has a radial field outside and none inside. 



so 



thus 



(1 - ^l')dPn/dn = !^^±ll(P„_i - P„+i) 



^ /a\2m+l X2n+1 p A ^ 

Notice that X2n+i/(4n + 3) is precisely our so 

E (f)'"^' ^2n+l(i^2n " P2„+2) r > a 
E i^Y"^' ^2„+l(P2n - P2n+2) T < a 



(38) 



A similar change of potential shows that the electric field 
can be written as V0 x VC where 



L = / r X V(t)V^.VCd^r = [ V^.VCd^r , 

Aire J VKc J 

where R x Wcj) = z and R = {x, y, 0) radial components 
of r X average to zero around a circle about the axis, 
by symmetry. 

Now V$.VC = div(CV$) - CV^$ and is zero 
both within and outside the sphere. Hence the angular 
momentum in the field within the sphere is 



U = £j div(cv$)dV = £ / cv$.d^^ , 



2c 



dr 



a d/j. , 



where the gradient d^/dr is evaluated on the inside of 
r = a. Similarly for the angular momentum external to 
the sphere 



Lo = + 



'2c J '^~dr 



a+ 



SO the total angular momentum depends on 
(fU-(f)a- = -4^^a^dwefind 



C = q 



, ^ 1 \ a'^uil + u) 1^ 2\ 



A 



Ar2 32 



+ 



[(r-iav/^)2]i/2 



r > a (39) 



The above expression is useful as the contours of constant 
C give the electric lines of force. Whereas these are radial 
for a non-rotating sphere they are more interesting for 
high rotation speeds and are plotted for u = 0.9, 7 = 19 
as Figure 1. 



L = Lo -h Lj = — z27ra^c 
Inserting the expressions for ( and a from the last section 



r + l 

-1 J Qadii . 



L = 



2c 



/: 



^2n+l(f'2n — -P2n+2) 

,n=0 

00 

X ^ (4m -I- l)$2m-P2mrfM 



m=0 
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Most of the terms vanish due to the orthogonality of the 
P„(/i) and we are left with (defining = 0) 



similarly there is no magnetic quadrupole moment but 
the electric quadrupole is 



But by equation (7) this is just (remembering $0 = 1) 



L = z— (E + w?7i) where I] = y^$ 
cuj ^-^ 



2 

2m 



(40) 



These formulae are exact without any approximation. In 
a rather similar manner one may derive exact formulae 
for the electical energy ee = J {E'^ /8TT)dPr and the total 
both inside and out gives 



1 f 

e. = -|^ c,4..2.a^,,= |_^$^„.|_(l + S) 

(41) 

For the magnetic energy we have 

= — / (curlA)2dV = — / div(A x B)dV 
= ^ JiAxBUS- J{AxB)odS 
= j ^^[Bo - Bj](,27ra2d/i 



asm( 



- A-Kauj sin 9d^ 



r + l 

/ ^(4to+ l)$2mf2md/^ 
J-1 ™ 



—Loq 

-— (E + cjryi) = 

2 a 2 a 



(42) 



Thus 



e = ee + e„i = 



- + + 



E is given by (27) and '7i by (24) in our first approxima- 
tion. 



Qe 



similarly the odd moments are missing from electric field 
and the even moments are missing from the magnetic 
field. We shall not detail the (2n + 1)*^ order moments 
but we note that as u — + c they only fall as [n + 3/4) ~" 
with a close to 0.2. This slow fall leads to a divergent 
total field energy as u — > c. This may be seen by com- 
paring the series in the sum E with that for Riemann's 
zeta function ^(1) which is infinite. 

The gyromagnetic ratio is fJ-m/L and the g- factor 



2mc fim 



so 



(^ + E + f^^i) 

^ (6/mc2) (ry, + j^T) 

This formula along with (|40|) assumes that all the angular 
momentum is in the moment of the Poynting vector of 
the electromagnetic field. When the system is highly 
relativistic ^ dominates so 
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2Vi 



9 / mc 



{e/mc^) 8 



where e is the energy in the field. Notice that the presence 
of non-electromagnetic mass will increase g above 9/8 
if we have only electromagnetic angular momentum. 



VII. OLD FASHIONED ELECTRON MODELS 

A century ago it would have been natural to apply the 
foregoing results to Mie's 'model' of the electron. Now ev- 
eryone knows better (and some like me understand less). 
However in the spirit of those times let us set the charge 
equal to e, the dipole moment equal to eh/2mc and the 



(43) total field angular momentum equal to ^h. Then 



eh 
2mc 



-ea rji 



(45) 



VI. THE DIPOLE AND QUADRUPOLE 
MOMENTS 



Directly from our expressions for $ and x ^ series we 
may read off the dipole and quadrupole moments of the 
field. There is no electric dipole and the magnetic dipole 
is 



-qaXi 



(44) 



Lc 



1 



2e2 ' 2ui/2 
and the total electromagnetic energy is then 




mc 



(46) 



(47) 
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From H46|l we see that fixing the angular momentum 
and charge determines the relativistic factor. Since Vi is 
never large and is only 0.375 in the extreme relativistic 
limit where S becomes large the fact that the left hand 
side is ~ ^ shows us that we need the highly relativistic 
regime with u close to 1 and S ~ 68. In that regime the 
terms in E dominate both numerator and denominator of 
the right hand side of (|47|l so e/mc^ ~ ^Vi = jg where we 
use the limiting value of Vi for this extremely relativistic 
case ^ 265, 7 ~ 530. Returning to eauation l|lS|) 
we find a = ^h/mc so the radius of the sphere is a little 
over ^ Compton wavelengths. The fact that at this 
large radius over half the rest mass is in the field energy 
shows that it was neglect of the effects of spin and charge 
polarization that led to the smaller radius /rac^ being 
considered classically. 

Is there anything to be learned from these naiive calcu- 
lations? Perhaps yes: Firstly polarization charges in the 
highly relativistic regime are much higher than the net 
charge. This helps classical motions to give angular mo- 
menta much higher than g^/c. Carter already pointed 
out that g — 2 for the Kerr metrics in relativity but 
did not draw attention to the polarization (Lynden-Bell 
2003). g = 2 is a generic property of all the conformasta- 
tionary metrics in General Relativity as well as the Kerr 
Metrics and Pfister & King (2002) suggest it is a rather 
general result. When the charges are uniform and fixed 
to the sphere the g factor varies from 3/2 to 11/6 as v 
increases (see Appendix I) . Secondly the value of the elec- 
tromagnetic angular momentum depends on / c and on 
how relativistically the system rotates but at given v/c; 



it does not depend on the radius of the configuration. 
The electromagnetic field energy and the dipole moment 
do depend on size, the latter gives a size related to the 
Compton wavelength but although this is two orders of 
magnitude larger than the classical radius of the electron 
the polarization due to rotation is so strong that more 
than half the rest mass can be electromagnetic even at 
this large radius. 

Finally the fact that tjrac? is quite near i suggests 
that a model with a somewhat different structure might 
yield exactly i . If the energy consists of an electromag- 
netic part e cx 1/a and a string with e' = Ta then equi- 
librium will occur when a d/ da{e -I- e') = i.e. e = e' in 
which case mc^ would be e -I- e' and e/mc^ = i. We note 
that a string loop will have no angular momentum about 
the axis of the loop so a purely magnetic angular momen- 
tum would then be correct. However my most recent disc 
models suggest that all the energy is electromagnetic. 
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APPENDIX A: THE RAPIDLY ROTATING 
UNIFORMLY CHARGED SPHERE 

The surface charge a = . ^ 2 is fixed in the surface. 
The resulting surface current is c'^^ilasmOa = J. We 
show that a field inside of 2fia~^ and an external field of 
a dipole (3/i..ff — /x)/r^ satisfies the boundary conditions 
on the sphere. 



sphere from inside to outside. 

Secondly [S^]™* = fia-^sm9{l + 2) this has to equal 
AttJ = c~^fta~^qsmO so the boundary condition is sat- 
isfied provided fi = ^qila^/c = ^qa^- The electric field 
is E = (/f/?-^ outside and zero inside. 

The angular momentum in the field is found from 
Poynting's vector 



i/rx(ExB)rfF.^/rx(^Ex-^)dy 



Attc J \ r"^ 

L = ^fiq J J{2cos'^9-3cos^9 + l)^r^sm9d9dr 
2 , 2 q^ V 



The energy in the electric field is ^g^/a = eg. The energy 



Firstly Br — 2^,a cos 9 and is continuous across the in the magnetic field is Cm where 
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3/ 3™ 



J (Scos^e + l) r-^dV 



^2 2 g2 ^2 
fl'^ 9 a 

The gyromagnetic 

ratio is n/L = ^ac/q. If all the angular momentum is in 
the field then the ^-factor is 



Thus e = e„ 



2 a 



2 

9 ^ 



2mc fi mc a 3 



Notice two results of this calculation 



2 f_ 

9 c2 



1. ^ ^ 



^ I for ?; ^ c so putting L = ^ and g = e there 
could not be an angular momentum in the field as 
great as |. The maximum seems to be short by 
the large factor '--^ | x 137. But the conducting 
sphere behaves differently. 

2. The g factor is always less than 2 if the energy 
is entirely electromagnetic but it is always greater 
than 3/2 and when v = c it becomes 11/6. Of 
course if the total energy is not all electromagnetic 
this result will change. If we fix the relativistic 
factor then e oc 1/a if the other energy e„ oc a" 
then the condition of equilibrium gives e = ne„ so 
mc^ = e + e„ = (1 + l/n)e. Then the factor ^ = 
(1 + 1/n). Cases could be made for n = 1 a string 
holding the electron together, n = 2 a tension over 
an area or n = 3 a tension throughout a volume so 
the factor might be 2, 3/2, 4/3 or of course 1 and it 
is of some interest to decide between these factors. 
However in the light of (1) above the insulating 
sphere with charge uniformly fixed to its surface is 
a bad model. 



APPENDIX B: 



SOLUTION BY GENERATING 
FUNCTIONS 



The equations to be solved are jSJ and ^ and we 
rewrite the latter in the form 



5Z(l+'5„)'72„+l = J$2. 



We multiply by C/™ and sum from m = to oo. We 
define the generating functions $(C/) and '?([/) by 



$([/) = ^$2mC/™ and %U) = 'WiC/" , 



m=0 



and so discover 



oo oo 



m— n—m 



Now 



and 



thus 



oo oo oo n 

EE-EE 

m— n—m n— m— 



Y u"' = (i-c/"+^)/(i-[/) , 

m=0 



^(1 - C/"+l)^2„+l + 5^(1 - f/"+l)^„'72„+l 



Ln=0 



= (1 - [/)-$([/) 



SO 



?7(1) - U7j{U) + ^(1 - [/"+i)(5„??2„+i 



^-(l-UMU) ■ (Bl) 

We now perform similar operations on equation (jsj 
multiplying by U"^ and summing from m = to oo. 
Reversing the order of summation 



Here we give a second method of solving the mathe- 
matical problem which allows a greater accuracy to be 
obtained by providing a first order perturbation treat- 
ment to the terms totally neglected in the first method. 
Also this second method provides the basis for solving 
the problem of the relativistically rotating disk which we 
treat elsewhere. We have refrained from replacing the 
first treatment by the second because the second, though 
more powerful, is less direct and guidance from the first gg 
method proved vital in understanding some of the steps 
essential to the development of the second. 



J2 E '^"^^n = LorjiU) 



n— 1 m—O 



J2 U"' = (1 - t/")/(i - U) 



m=0 



$(1) - $([/) = -uj{l - U)ri{U) 



(B2) 
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Notice that the lack of an n = term in the sum does 
not matter as it cancels out between $(1) and $(?7). 
Putting U = in the result gives 

$(1) - 1 = -c^ 771 , (B3) 
so we get an expression for $(?7) in terms of rj 

$(f/) = uj{l - U)'n{U) + 1-Lur]i . 
Inserting this into HB1() and solving for ri{U) 



rjiU) 



'7(1) + (1 - C/) ( - ^ 



^(l-C/"+i),5A.+i 



(B4) 



where Q is the quadratic obtained from H12|l by writing 

t = -U 



Q = U^ + 2{2uj-^ ~1)U+1 = {U + u){U + it-i) . 

In our former approximation we neglected the (5„ for n ^ 
2 and kept only and 5i. Doing that in l|B4(l leads, via 
the method we shall shortly explain, back to our former 
solution which gave ^2n+i = C(— m)" for n ^ 2. We shall 
adopt this approximate form of solution only for the small 
terms with <5„, 7i ^ 2 in the sum in (|B4|I but we shall 
not assume that the value of C is exactly that given in 
equations (|25|l and (|26|l and we shall derive refined values 
for the '?2n+i- Thus the sum in (|B4p becomes 



^(1 - C/"+i)(5„'?2„+i = (1 - U)5a^i + 



oo 

+ (1 - C/2)<5ir/3 + C^(l - U^+^)6n{-uY (B5) 

2 

We now define a function W{t) 



oo 1 1 



Writing tip' for t we find 



and 



d_ 

dw 



w'^ — {wW) 
aw 



w 
1 



w 



4(1 + w2) 



d/dw{wW) = T-^ln(l + w'^) 

and integration by parts then gives, again keeping W 
finite at u; = 



1 1 

^ 4w Jo 



dw 




Replacing w'^ by t and taking account of the possibility 
that t may be negative we find the alternative forms of 
expression of a single analytic function 



Wit) 



Notice that the series expansion of the final functions 
only involve integer powers of t so the square roots 
are only involved in the formal expression in terms of 
tan~^&ln. W{t) is analytic in the complex t plane in 
\t\ < 1. Our former approximation retains only So and i5i 
so replaces W by the linear approximation Wl = 
At i = 0, W = Wl = 1, aw = 
At t ^ 1, W ^ ^ - iln2 = 0.1097, Wl 
0.1042, AW = 0.0055 



48'"- 



W 



i In 2 = 0.1733, Wl = 0.1458, AW = 



At t = -1, 
0.0275 

So while that approximation for W remains good to 5% 
in the more physical range of positive t, the accuracy falls 
to only 13% at t = —1 and our arguments will involve 
W there. It is therefore of importance to get the more 
accurate result to which this Appendix is devoted. Re- 
expressing (|B5ll in terms of W 



^(1 - C/"+i)<5„?72„+i = (1 - U)6oiVi -C) + 



+ (1 - U^)6i{V3 + Cu) + C[W{u) - UWiuU)] 
so our expression (|B4|) for ri{U) becomes 



{77(1) + (1 - C/) 



Integrating and ensuring convergence at w = gives us analytic i.e. within \U\ 



(1 - C/2)(5i(f?3 + Cu) + C[W{u) ~ UW{uU)]} .(B7) 

'72i+i is the coefficient of C/' in this re-expressed as a 
power series in U. However in our solution pi(l of the 
recurrence relation (jlO|l we had to choose the convergent 
solution and throw away the divergent one and only then 
did we have enough equations to yield a solution. We now 
have to find the analogous step in our new approach. The 
numerator of r](U) is analytic everywhere that W{uU) is 



Thus any poles that 
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can give us divergences in the high order coefficients of 
[/' can only come from the zeros of Q. One such pole 
\s &t U — —u~^ which is outside the \U\ — 1 circle, 
inside which rjiU) should be analytic, but the other is 
within that circle at [/ = u. It is this pole that leads to 
the divergent solution of the recurrence relation whose 
coefficient D had to be chosen as zero in Hll|) . Thus 
the requirement of convergence in our coefficients ^2n+i 
can be replaced by the requirement that r?(C/) should be 
analytic with no poles within the circle \U\ = 1. This in 
turn translates into the requirement that the numerator 
of rj{U) must vanish aX U = —u so 



7/(l) + (l + u) 



+ '5o(^i-C) 



+ (1 - u^)5i{Tls. + Cu) + C[W{u) + uW{-u^)] ^ 

without first having the solution by the recurrence rela- 
tionship it would have been difficult to see this require- 
ment. We now use this form of 77(1) in (|B7p . At the same 
time we replace i^'^ ,Q by their values in terms of u given 
in equation (|13|l . Then 



from the first term which gives the asymptotic expan- 
sion CS (—[/«)" as expected from the asymptotic form 
of the recurrence relation. This justifies the nomencla- 
ture of calling its coefficient C which was shorthand for 
(cf llBlOll 'l 



C = Vi + [l + ufu-^ 



Si{V3 + Cu) + CU 



W{-u^) - W{-1) 



1 



(B12) 



Now from the definition of ri{U) we have ??'(0) = V3. 
Applying this to (|B11|) we find the third equation which 
together with ljB9|l and ljB12|l we solve for Vi, ^3 and C 



V3 = -Cu + Cil 
Wi-u^) ~W{-1) W{-u^) ~WiO) 



u? X 



-Cu + Cu{l + ufgo (B13) 



1 



The resulting value of '73 -I- Cu we substitute into (|B12p 
to obtain 



r]{U) = (l + Uu) 1 -uryi + iM^(l + 

+ (1 + uf {-So{Vi ~C) + {u- U)Si{V3 + Cu)- 
- C[uW{-u^) + UW{uU)]/{u + U)}] (B8) 

Now 77(0) = Vi hence 

Vi = -uT]i + ^u^/^{l + u) + {l + uf{-So{Vi~C) + 
+ u5i{V3 + Cu)~CWi-u^)} (B9) 
We use this expression to simplify IIB8I) into 



r^{U) = {l + Uu)-^{Vi + {l + ufUx 

'^(-^^)-^("^)n(Bio) 



-Si{V3 + Cu)+C- 



U + u 



Notice that the term involving the Ws has no singular- 
ity when U = —u and that term in the overall numerator 
takes the value 

-ulWi-u^) - W{-1)]/{1 - u^) 

when U = —l/u. This is finite for u < 1. Thus isolating 
the term with a pole at J7 = — 1/u we have 



c 



(1 + uf 



1 + Uu 
U 



Wi-u^)-W{uU) 

U+u 



5i{r]^ + Cu) + C{l + uf X 

W(-u'^)-W(-l)~ 



l + Uu 



(Bll) 



where the numerator of the last term is zero at [/ = — 1/m 
so it has no pole there either. The first pole comes 



Vi = C{1 - (1 



(B14) 



where 51 = — — ^-M^i — 9o = 9i 
when only 5q and 5i are retained gi = —5i and go is zero. 

We now use (|B13|) and (|B14p in (|B9|1 to obtain an 
expression for C 



W{-u^)-W(0) 



c 



,1/2 



/A 



where 



A 



= l-{l + uf{gi[l + {l + u)5o] + 

+ {l + u)5igo[l + u + y?) + {l + uf5o]} + 

+ {l + u)W{~u'^) (B15) 



which reduces to H26|) in the appropriate limit. With 
"Hi and C now determined 'ri{U) is known from IjBlip . 

We are now in a position to obtain a more accurate 
value for the surface density of charge at the poles and the 
equatorial velocity at which it changes sign. According 
to ^ we need 



E 





(4n + 1)$2„ - 4$'(1) + $(1) 



But by (in $'(1) = -ujr]{l) 
and by ijlHIl $(1) = I - LU-qi 
so ap ^ - u;[V, + Ml)]} 

Evaluating this expression as a function of u we find 



that CTb 



when u = 0.455 which is in better than 



expected agreement with the crude calculation of section 
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4 which gave u = .454 or 93% of c for the equatorial 
speed. 

We now have our sokition: C is given by 1B14|I and 
(|B15p . Let us recap the notation uj = ila/c — ^ 



In practice the most interesting region is the relativistic 
region, u = ^ gives 7 = 3, u = 0.8 gives 7 = 9 while 
u — 0.9 gives 7 = 19. Evaluating W{—u'^) for such values 
we find 



= 2cj" 



1 - 



1 
2 

-1/2 



1 



so 7 



W{t) is the function defined by IjBSp 
g\ is the expression u^) — 1)]/(1 

30 is gx+u--'\W{-v?)-Wm 
70 = W^(0) = i , <Si = ^ 

Thus all terms in C are known once w or u are given. 
(|BT3|) then gives '?i and lj5T2)l gives n^. 
The function 77(C/) is now given by ljB7|) with all terms 
therein known. 

The coefficient of [/" in the expansion of ?7(?7) as a 
power series in U is our refined solution for V2n+i but 
with the first two Vi and V3 now explicitly known the 
others may also be found by direct repeated application 
of the exact recurrence relation 110|l e.g. 



V5^C{u^ + (l + u)^[g,+Si-go{l + u)^]} 



Wi-u"^) = — [(1 - u) ln(l -u) + {l+u) ln(l + u)] 



8u 



Near u — 1 gi 



W^(-l) = iln2 = .173, iy(0) = i 



^8u2(l+«) ■ 



x[ln(l -u) - (ln2- 1) + (1 - u)(ln 2 - i)] 
iu"'^[(l - M)ln(l - u) + (l + M)ln(l + u)-u2] 



the dominant term is the ln(l — u) term in gi. 

Notice that near u = 1 these terms become large which 
invalidates our perturbation theory beyond 7 ^ 20. The 
7 > 20 region is treated differently in the main body of 
the paper. 



